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Abstract: The time-harmonic electromagnetic responses of (a) a bilayer made of an 
epsilon-negative layer and a mu-negative layer, and (b) a single layer of a negative phase- 
velocity material are compared. Provided all layers are electrically thin, a restricted 
equivalence between (a) and (b) exists. The restricted equivalence depends on the lin- 
ear polarization state and the transverse wavenumber. Implications for perfect lenses and 
parallel-plate waveguides are considered. 
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1 Introduction 



This communication is inspired by the ongoing spate of papers published on the inap- 
propriately designated left-handed materials which are macroscopically homogeneous and 
display negative phase velocities in relation to the time-averaged Poynting vector, but 
are not chiral 1 . Nominally, such a material is deemed to possess a relative permittivity 
scalar e r = e' r + ie" and a relative permeability scalar \i r = fi' r + in", both dependent on 
the angular frequency oj, such that both e' r < and \J r < in some spectral regime. 2 
Originally conceived by Veselago Ves, these materials have been artificially realized quite 
recently SSS,PGXKT. Their fascinating electromagnetic properties can have technological 
implications of massive proportions Pen, but those implications remain speculative at this 
time. Although these materials have been variously named LMW, we prefer the name 
negative phase-velocity (NPV) materials as the least ambiguous of all extant names. 

Using transmission-line analysis and lumped parameters, Alu and Engheta Engh have 
recently suggested a new route to realizing NPV materials: Take a thin layer of an epsilon- 
negative (EN) material: it has a negative real permittivity scalar but a positive real 
permeability scalar. Stick it to a thin layer of a mu-negative (MN) material, which has 
a negative real permeability scalar and a positive real permittivity scalar. Provided the 
two layers are sufficiently thin, the paired EN-MN layers could function effectively as a 
NPV material. The clear attraction of this scheme is that EN and MN layers are easier 
to manufacture, very likely, than the NPV materials fabricated thus far SSS,PGXKT. 

Our objective here is to examine the suggested scheme using continuum field theory, and 

to establish a restricted equivalence of an EN-MN bilayer to a NPV material. Implications 
2 The condition for the phase velocity and the time-averaged Poynting vector to be oppositely directed 
is (|e r | — e.' T ) (|Mr| — aO > which permits — more generally — e' r and/or y! r to be negative MLW. 

An exp(— iuit) time-dependence having been assumed here, e" > and /i" > at all ui > for all passive 
materials. 
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for parallel-plate waveguides Engh and perfect lenses P2 are deduced therefrom. 

A note about notation: Vectors are in boldface, column vectors are boldface and en- 
closed in square brackets, while matrixes are denoted by Gothic letters; eo and are the 
free-space permittivity and permeability, respectively; and k = w(eo^o) 1 ^ 2 is the free- 
space wavenumber. A cartesian coordinate system is used, with u^,, u y and u z as the unit 
vectors. 



2 Bilayer Theory in Brief 



Consider the layers < z < d a and d a < z < d a + dj,. Their constitutive relations are as 
follows: 

D(r) =e e a E(r) 



B(r) = |i |i a H(r) 



< z < d a , 



(1) 



d a < z < d a + db ■ 



(2) 



D(r) =e e b B(r) 
B(r) = w t H(r) 

The constitutive parameters present in the foregoing equations are complex-valued with 
positive imaginary parts (as befits any passive medium). The two half-spaces on either 
sides of the bilayer are vacuous. 

Without loss of generality, the electromagnetic field phasors everywhere can be written 
as Kl 

E(r) = e(z)e" 



> , — oo < z < oo . 



(3) 



H(r) = h(z) e lKX 

where the transverse wavenumber k € [0, oo). The fields inside the bilayer must follow 
the 4x4 matrix ordinary differential equation Kl^Ll 



^ [f (Z)] = iV(z) [f (Z)] , 0<Z<d a +d b . 



(4) 
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In this equation, [f(z)] = col e x (z), e y (z), h x (z), h y (z) is a column vector, while the 4x4 
matrix function Viz) is piecewise uniform as 



V{z) 



where 









V a , 

n, 






< z < d a 
d a < z <d a + d b 





(5) 













(6) 



The only nonzero elements of the matrixes V a ,b appear on their antidiagonals, of which 
the (2,3) and the (3,2) elements are relevant to s-polarized fields, and the (1,4) and the 
(4, 1) elements to the p-polarized fields. 

The solution of diffeq is straightforward, because the matrix V[z) is piecewise uniform 
Hqch. Thus, the algebraic relation 



[f (d b + d a )} = e iVbd » e iVada [f (0)] 



(7) 



is sufficient to solve both reflection/transmission problems as well as guided-wave prop- 
agation problems. The two matrix exponentials on the right side of sys cannot be inter- 
changed — unless the matrixes V a and V b also commute, which is possible with dissimilar 
materials only in quite special circumstances L1 4 L2. 



3 Analysis 



Matrixes V a ,b have ±y ^e a ^ a f) — k 2 = ±a ayb as their eigenvalues. Provided that \a a)b \ d a , b <C 
1 (i.e., both layers are electrically thin), the approximations 



€? V «*< b ~ l + iV a ,bda,b (8) 



can be made, with X as the 4x4 identity matrix. Then 

e iV h d h e i,V a d a _ J- + iVada + _ e i,Vad a e iV b d b ^ (g) 

and the two layers in the bilayer can be interchanged without significant effect Reese. 

Let us now consider a single layer of relative permittivity e eq , relative permeability fi eq 
and thickness d eq . Quantities V eq and a eq can be defined in analogy to V a and a a . Two 
thickness ratios are defined as 

Pa, b =^>0, (10) 

in order to compare the single layer with the previously described bilayer. There is no 
hidden restriction on the non-negative real numbers p a and p\>. 

Provided that \a eq \d eq <C 1 (i.e., the single layer is electrically thin as well), the ap- 
proximation 

e iVeqdeq -1 + iVeqdeq (11) 

can be made. Equations ab and eq permit us to establish the following equivalences 
between a bilayer and a single layer: 

(i) s-polarization: The only nonzero field components are E y , H x and H z . Therefore, 
the equality of the (2, 3) elements of p a V a + PbPb and V eq , and likewise of the (3, 2) 
elements, has to be guaranteed for equivalence; thus, the equations 

PaVa + PbUb = Veq , (12) 

P a e a +Pbtb - (j-) (— + —') =^eq~ (-r-) — (13) 

have to solved for e eq and fj, eq . We conclude therefrom that, for a given value of k 
and subject to the thickness restrictions |o; a 5 e „| d a ^^ eq ~ 0.1, the bilayer and a single 
layer are equivalent with respect to the transformation of the x- and y- components 
of the fields from one exterior face to the other exterior face if 

e e q = Pa^a + Pb^b 
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PaPbiVa ~ Vbf iPa + Pb + l)(Po + Pb ~ 1) 



_HaVb(PaHa + PbUb) 
Pa^a + PbPb ■ 



PaPa + PbfJ-b 



(14) 
(15) 



(ii) p-polarization: The only nonzero field components being H y , E x and the equality 
of the (1,4) elements of p a V a +PbPb and V eq suffices, along with the equality of the 
(4, 1) elements of the two matrixes. For a given value of k and subject to the thickness 
restrictions |a a ,fe,eg| d aj b,eq ~ 0.1, the bilayer and a single layer are equivalent if 



^eq = Pa^a + PbUb 
K 



PaPbi^a ~ e b f . (Pa + Pb + l)(p a + p b - 1) 



teq — 



e a eb(p a ^a + Pbtb) 

Pa^a + Pbtb • 



+ 



Pata + Pbtb 



(16) 

(17) 



Clearly, the constitutive parameters of the equivalent layer are functions of both p a and p b ; 
and we must point out that the sum of these two ratios need not equal unity. Furthermore, 
except for normal incidence (i.e., n = 0), the constitutive parameters of the equivalent layer 
depend on the incident linear polarization state. Finally, the constitutive parameters of 
the equivalent layer change with the transverse wavenumber k. 

The foregoing equations can be manipulated to yield negative values of both e' and 



fj,' eq for either 



an EN-MN bilayer {e' a < 0, e' b > 0, fi' a > 0, n' b < 0} or 
a MN-EN bilayer {e' a > 0, e' b < 0, fjf a < 0, n' b > 0}. 



An EN-MN bilayer is equivalent to a NPV layer for both polarization states when k = 0, 
provided the condition 



Wb\ ^ Pa ^ < 



Ha 



> — > 



(18) 



Pb \e a \ 

holds true. The inequality deltal is applicable for a MN-EN bilayer also, if the subscripts 
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a and b are interchanged therein, i.e., 




(19) 



A further specialization of p a + = 1 leads to the inequality 



> Pa 



> 




(20) 



for EN-MN bilayers, and 



>Pb> 




(21) 



fi'b + Wa 



for MN-EN bilayers. 



The inequalities deltal~delta4 should be adequate for both s- and p-polarization states 



when k/J;q < 1. In general, however, a given EN-MN (or MN-EN) bilayer is equivalent 
to a different NPV material for a different linear polarization state and/or transverse 
wavenumber. Thus, the equivalence between an EN-MN (or a MN-EN) bilayer and a 
NPV layer is restricted. 

The restricted equivalence has an interesting implication for perfect lenses P2. A perfect 
lens of thickness d > is defined by the fulfillment of the condition [f(d)] = [f(0)] for 
all co and k. Because of dispersion and dissipation, at best, this condition is fulfilled 
approximately. Let us imagine that the condition is fulfilled by some NPV constitutive 
parameters for some to and all \k\ < k. Then, the implementation of the acceptably 
imperfect lens as a cascade of thin EN-MN (or MN-EN) bilayers would require that the 
successive bilayers have different constitutive parameters and that the entry as well as the 
exit faces be curved, and even those steps may not suffice. In contrast, k is fixed for any 
single-mode parallel-plate waveguide, and so is the range of operating frequencies; and 
the emulation of a NPV material via EN-MN (or MN-EN) bilayers may not be onerous. 



7 



References 

[1] Lakhtakia A, McCall MW, Weiglhofer WS: Brief overview of recent developments on 
negative phase-velocity mediums (alias left-handed materials). Arch. Elektr. Uber. 
56 (2002) 407-410 

[2] McCall MW, Lakhtakia A, Weiglhofer WS: The negative index of refraction demys- 
tified. Eur. J. Phys. 23 (2002) 353-359 

[3] Veselago VS: The electrodynamics of substances with simultaneously negative values 
of e and fi. Sov. Phys. Usp. 10 (1968) 509-514 

[4] Shelby RA, Smith DR, Schultz S: Experimental verification of a negative index of 
refraction. Science 292 (2001) 77-79 

[5] Parazzoli CG, Greegor RB, Li K, Koltenbah BEC, Tanielian M: Experimental veri- 
fication and simulation of negative index of refraction using Snell's law. Phys. Rev. 
Lett. 90 (2003) 107401 

[6] Pendry P: Electromagnetic materials enter the negative age. Physics World 14 (9) 
(2001) 47-51, September issue 

[7] Alu A, Engheta N: Mono-modal waveguides filled with parallel epsilon-negative and 
mu-negative paired metamaterial layers. 2003 IEEE MTT-S Int. Microw. Symp., 
Philadelphia, PA, USA, June 8-13, 2003 

[8] Pendry J: Negative refraction makes a perfect lens. Phys. Rev. Lett. 85 (2001) 3966- 
3969 

[9] Krowne CM: Fourier transformed matrix method of finding propagation character- 
istics of complex anisotropic layered media. IEEE Trans. Microw. Theory Tech. 32 
(1984) 1617-1625 



8 



[10] Lakhtakia A: On perfect lenses and nihility. Int. J. Infrared Millim. Waves 23 (2002) 
339-343 

[11] Lakhtakia A, Sherwin JA: Orthorhombic materials and perfect lenses. Int. J. Infrared 
Millim. Waves 24 (2003) 19-23 

[12] Hochstadt H: Differential Equations — A Modern Approach, Dover Press, New York, 
NY, USA 1975 

[13] Reese PS, Lakhtakia A: Low-frequency electromagnetic properties of an alternating 
stack of thin uniaxial dielectric laminae and uniaxial magnetic laminae. Zeit. Natur- 
forsch. A 46 (1991) 384-388 



9 



